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Abstract We study the electronic transport in quasiperiodic separable tight-binding 
models in one, two, and three dimensions. First, we investigate a one-dimensional 
quasiperiodic chain, in which the atoms are coupled by weak and strong bonds 
aligned according to the Fibonacci chain. The associated tZ-dimensional quasiperi- 
odic tilings are constructed from the product of d such chains, which yields either 
the square/cubic Fibonacci tiling or the labyrinth tiling. We study the scaling be- 
havior of the mean square displacement and the return probability of wave packets 
with respect to time. We also discuss results of renormalization group approaches 
and lower bounds for the scaling exponent of the width of the wave packet. 



1 Intoduction 

Understanding the relations between the atomic structure and the physical properties 
of materials remains one of the elementary questions of condensed-matter physics 
further emphasized by the discovery of quasicrystals 1 1 1. Quasicrystals are charac- 
terized by a perfect long range order without having a three-dimensional transla- 
tional periodicity. The former is manifested by the occurrence of sharp spots in the 
diffraction pattern and the latter in the occurrence of rotational symmetries forbid- 
den for conventional crystals. Already in the 1970s, works by Penrose and Ammann 
showed that the Euclidean space can be filled gapless and non-overlapping by two 
or more tiles which are arranged in a nonperiodic way according to matching rules. 
It turned out that these tilings are suitable to describe the structure of quasicrystals. 
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Experimental studies revealed rather exotic physical properties. For instance, 
quasicrystalline surfaces are anti-adhesive in combination with a high level of hard- 
ness making them suitable for the production of coatings for medical equipment, 
engines, cookware, etc. Further, they possess a low thermal and electrical conduc- 
tance although they contain a high amount of well-conducting elements. For one- 
dimensional quasicrystals many numerical studies helped to establish a better un- 
derstanding of the physical properties [2, 3, 4|. Further, several exact theoretical 
results are known in one dimension today |5 |. However, the characteristics in two 
or three dimensions have been clarified to a much lesser degree because numerical 
studies are usually restricted to small approximants. 

To address this problem we study models of c/-dimensional quasicrystals with a 
separable Hamiltonian in a tight-binding approach |6|. This method is based on the 
Fibonacci sequence, which describes the weak and strong couplings of atoms in a 
quasiperiodic chain. After a iterations of the inflation rule ^ = {s — > w, w — > ws}, 
we obtain the ath order approximant of the Fibonacci chain. The length of an 
approximant '^a is given by the recursive rule fa = fa-i + fa-2 with /o = /i = 
1 . Further, the ratio of the lengths of two successive approximants approaches the 
golden mean T for a — > °o. Solving the time-independent Schrodinger equation 

H|f') ^t,_u1'/-i +tij+i1'/+i =EX (1) 

for the Fibonacci chain, we obtain discrete energy values E' and wave functions 
jf') ~ Lflt' "^Z 10 represented in the orthogonal basis states |/) associated to a ver- 
tex /. The hopping strength t in the Schrodinger equation is given by the Fibonacci 
sequence with = s for a strong bond and tw = w for a weak bond (0 < w < s). 

The (i-dimensional separable quasiperiodic tilings are then constructed from the 
product of d quasiperiodic chains which are perpendicular to each other. For this 
setup two special cases are known for which the systems becomes separable |l6l: 

• Hypercubic Tiling J^/^ — This tiUng corresponds to the usual Euclidian prod- 
uct of d linear quasiperiodic chains, i.e., only vertices connected by vertical and 
horizontal bonds interact as shown in Figure [T] The electronic structure and the 
transport properties are understood quite well for these systems EH. 

• Labyrinth Tiling ^^'^ — Only coupling terms to neighbors along the diagonal 
bonds are considered (cf. Figure [TJ, where the bond strengths of this tiling equal 
the products of the corresponding bond strengths of the one-dimensional chains. 

The eigenstates of these tilings in d dimensions can be constructed from the eigen- 
states of d one-dimensional chains. In both cases the wave functions of the higher- 
dimensional tiling are constructed as the products of the one-dimensional wave 
functions, i.e., = ^fj^^^L °^ I'l^'^'m ■ ■ ■I'f ■ The superscripts s = {i,j,...,k) 
enumerate the eigenvalues E, and r = {l,m,...,n) denotes the coordinates of the 
vertices in the tiling. The energies are treated in a different way. In d dimensions 
they are given by E^ = E'--' =E^' +E~^ + . . . +E^'^ for the hypercubic tiling and 
by E^ = E''j'---'^ = E^'E^-i . . .E'^'^ for the labyrinth tiling. This approach allows us to 
study very large systems in higher dimensions with up to lO"^ sites. 
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Fig. 1 Two-dimensional 
square tiling J^^d (jg^j^ 

and labyrinth tiling 
(right) constructed from two 
Fibonacci chains . 

2 Quantum Diffusion 

To obtain a deeper understanding of the connections of the electronic transport 
and the quasiperiodic structure of a system, we study the time evolution of wave 
packets |r(ro,f)) = ErejSf ^('"OiO |r) initially localized at a position Tq which are 
constructed from the solutions ]^(ro,f) = Ls^ro^r^ °f time-dependent 
Schrodinger equation. The temporal autocorrelation function of a wave packet 
equals the integrated probability to be at the position ro up to time f > E), i.e., 
C(ro,f) = J Jq \YrQ{ro,t')\^dt' . We denote the integrand as return probability 

P(ro,f) = |rr„(ro,f)|'. (2) 

Another quantity often considered for the description of the electronic transport 
properties is the mean square displacement of the wave packet (also called width) 

fl'(ro,f) = ^ErG^|r-ronrr(ro,f)|2. (3) 

The wave-packet dynamics reveal anomalous diffusion with c/(ro,f) °^ f/^(''o), and 
the electronic transport is governed by the wave-packet dynamics averaged over 
different initial positions ro, i.e., d{t) = {d{ro,t)) °^ . The exponent j3 is related to 
the conductivity a via the generalized Drude formula, where ji —Q corresponds to 
no diffusion, j3 = 1/2 to classical diffusion, and j3 = 1 to ballistic spreading. 

The autocorrelation function is expected to decay with C(ro,f) °^ f^^t*"") Jgl. The 
exponent 5(ro) is equivalent to the correlation dimension of the local density 
of states (LDOS) p{r(),E) of a system |8|. In one dimension 5 — ?> 1 corresponds 
to the ballistic motion of an electron. More information about the transport can be 
obtained by studying the return probability which shows a power-law behavior ac- 
cording to P(t) t^^ . The integration is only used to smooth the results. However, 
this leads to some disadvantages: For 5' ~ 1 one obtains an additional logarithmic 
contribution C(f) °^ ln(f)/f leading to 8 significantly lower than 1 |3|. Further, in 
higher dimensions 5' > 1 is possible for large coupling parameters w. This leads to 
the convergence of the integral in C(f) as shown in Figures|2]and[3][iJ|. 

Some typical results for the mean square displacement d{t) and the temporal 
autocorrelation function C(f) are shown in Figure |2] We average the results over 
different initial positions of the wave packet, i.e., we study the scaling behavior of 
d{t) = {d{rQ,t)) oc t-P and C(0 = (C(ro,0> °= f-^. The power-law behavior can be 
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Fig. 2 Averaged mean square displacement d{t) for the Fibonacci chain "^19 (left), averaged tem- 
poral autocorrelation function C(t) for "^20 (center) and the 2D labyrinth tiling (right), i = 1. 



observed over several orders of magnitude in time, before d{t) and C(f) approach 
a constant due to finite size effects. The corresponding scaling exponents j3 and 5 
obtained by least squares fits are compiled in Figurejs] We find 5' > 1 for w > 0.55s 
in two dimensions and w > QA3s in three dimensions. Consequently, the scaling 
exponent 5 slowly approaches 1 for this threshold w as indicated by the dashed 
lines [9J. This is in reasonable agreement with the known transition to an absolute 
continuous energy spectrum for w^j^ « 0.6s and ss 0.46s QITOl. The exponent /3 
indicates anomalous transport for the hypercubic and the labyrinth tiling for all w < 
s. Only for w — s ballistic transport is found. Further, we observe that the exponent 
is about twice as large as the one-dimensional exponent 5[^. However, the results 
for the three-dimensional labyrinth tiling do not entirely fit into this scheme. While 
for w < 0.6s the expression S^^/d is only slightly smaller than the exponent 5[^, 
the behavior changes completely for large values of w. In the latter case the scaling 
exponent 5i^^ becomes almost constant and approaches 5i^^{w = s) « 1.75. 

For the hypercubic tiling the relation ~ ^mI^ holds, which also characterizes 
the transition to an absolutely continuous energy spectrum according to d8\i{w'l^) — 
1 [3|. While the numerical results for the 2D labyrinth tiling seem to satisfy these 
relations as well, there are significant differences in the wave-packet dynamics in 




coupling parameter w coupling parameter w coupling parameter w 



Fig. 3 Scaling exponents ji (left), 5 including the limit behavior (lines) for infinite systems (cen- 
ter), and 8' /d (right) for the Fibonacci chain and the labyrinth tilings in 2D and 3D for i = 1. RG 
results are shown by lines in left and right panel. 
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three dimensions, which are probably related to the different grid structure. In three 
dimensions for w s the labyrinth tiling approaches a body centered cubic lattice 
while the cubic tiling approaches a simple cubic lattice. 



3 RG Approach and Lower Bound for the ScaUng Exponent /3 

For the hypercubic tiling the scaling exponents of the width fulfill the relation jS^/d = 
/3id due to the separability of the time-evolution operator Numerical results suggest 
that this relation is also valid for the labyrinth. For weak couplings (w <C s) it is 
possible to find a connection between the quasiperiodic structure of the Fibonacci 
chain ifTTI or the labyrinth tiling lfT2l [T3]| with the wave-packet dynamics by an 
RG approach proposed by Niu and Nori 1141 . The results show that the scaling 
exponents j3 for different dimensions approach each other for w — ^ but are not 
identical lfT2l [131 (cf. Figure [3]). The hierarchic structure of the RG enables us to 
describe the properties of the wave functions and energy spectrum of the Fibonacci 
chain for w <C i fl4^'T5'|. This can be used to derive analytical relations, too, for the 
mean-square displacement 1 15 1 and the return probability (cf. Figure[3]l lfT3l[T6l . 

The computation of the scaling exponent j3 is computationally very expensive. 
However, one can make use of two general lower bounds for j3. As a rule of thumb, 
the wave-packet propagation is faster for smoother spectral measures. This was 
proved by Guarneri et al. with the lower bound j3 > jd based on the informa- 
tion dimension of the spectral measure of the LDOS |17|. Ketzmerick et al. 
showed that the spreading of a wave packet in a space of reduced dimension 
(i.e. the correlation dimension of the wave functions) is described by the bound 
j3 > D2/D2 as long as < 1, i.e., for singular continuous energy spectra lilSl . 
This result is based on the normalization condition and the known decay of the cen- 
ter of the wave packet according to f 2 . A comparison with the numerical results 
in Figure [4] shows that all inequalities are clearly satisfied. However, in the regime 
of an absolutely continuous energy spectrum, the spectral dimensions fulfill Dq = 1, 
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Fig. 4 Comparison of j8 with the lower bounds /d, /D2, and 8' /D2 for the Fibonacci chain 
(left) and the corresponding labyrinth tilings in two (center) and three dimensions (right) for .s = 1 . 
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and the two inequalities are no longer good lower bounds. We like to point out that 
the decay of the center of the wave packet is only described by the correlation di- 
mension D2 of the LDOS until the integral in C(f ) converges |3 |, i.e., as long as the 
exponent 5', which describes the decay of the center of the wave packet, is smaller 
than 1. Hence, for li-dimensional systems the bound by Ketzmerick et al. should 
be replaced by j3 > 8' /D2 in the absolutely continuous regime. For the hypercu- 
bic tiling this bound is as good as for the Fibonacci chain. This bound should hold 
also in general because it solely makes use of the normalization condition of the 
wave packet. We have checked in Figure]?] whether this relation is satisfied for the 
labyrinth tiling and found that this is a significantly better lower bound for w > Wth- 



4 Conclusion 



We have studied the wave-packet dynamics for separable quasiperiodic tilings. In 
higher dimensions 5^-1 indicates the existence of an absolutely continuous part in 
the energy spectrum rather than the occurrence of ballistic transport. To understand 
the diffusive properties of a system it is also necessary to compute the width d{t) 
of the wave packet and the return probability P{t) because in three dimensions the 
scaling exponents j3 hardly differ for the considered models although the scaling 
exponents 5' are rather different for large coupling parameters w. We also found 
good agreement between the scaling exponents 5' and j3 with the analytical expres- 
sions derived by an RG approach in the regime of strong quasiperiodic modulation. 
Further, the exponent 5' can be used to define a better lower bound for the scaling 
exponent /3 of the wave-packet width in the absolute continuous regime. 
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